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ABSTRACT 

The spectral distributions of Be/X-ray binaries in the Large Magellanic Cloud and 
Galaxy have been shown to differ significantly from the distribution of isolated Be 
stars in the Galaxy. Population synthesis models can explain this difference in spec- 
tral distributions through substantial angular momentum loss from the binary system. 
In this work we explore the spectral distribution of Bc/X-ray binaries in the Small 
Magellanic Cloud (SMC) using high signal-to-noise spectroscopy of a sample of 37 
optical counterparts to known X-ray pulsars. Our results show that the spectral dis- 
tribution of Be/X-ray binaries in the SMC is consistent with that of the Galaxy, despite 
the lower metallicity environment of the SMC. This may indicate that, although the 
metallicity of the SMC is conducive to the formation of a large number of HMXBs, the 
spectral distribution of these systems is likely to be most strongly influenced by angu- 
lar momentum losses during binary evolution, which are not particularly dependent 
on the local metallicity. 

Key words: stars: binaries - emission-line, Be: Magellanic Clouds. 
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1 INTRODUCTION 

Be/X-ray binaries are stellar systems in which a neutron 
star accretes material from a massive, early-type star with a 
circumstellar disc. Typically they have wide, eccentric orbits 
with X-ray outbursts occurring when the neutron star inter- 
acts with the circumstellar material of the Be star. Two pre- 
dominant types of outburst behaviour are observed: Type I 
outbursts of luminosities in the range 

10 36_ 10 37 ergg -i lagt 

a few days and generally occur close to periastron, while 
Type II outbursts last much longer, reach higher luminosi- 
ties (> 10 37 ergs _1 ) and show no correlation with orbital 
phase. 

The Small Magellanic Cloud (SMC) is home to an 
unusually large population of high mass X-ray binaries 
(HMXBs). Based on the relative masses of the Milky Way 
and the SMC, there are a factor 50 more in the SMC than 
one would expect. (Although the sample of HMXBs in the 
Milky Way is by no means complete, as can be evidenced 
by t he recent discoverie s of numerous obscured X-ray bina- 
ries l| Walter et al.l l2006). we do not expect the population 
to incre ase by a fact or of 50). Population synthesis simula- 
tions by iDravl (|2006l ) show that the reduced metallicity of 
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the SMC (~one-fifth solar) can lead to an increase in the 
HMXB population by a factor of three. The metallicity is 
particularly relevant in influencing the evolution of massive 
stars, due to the impact it has on stellar winds. The line- 
driven stellar winds in mas sive stars become weakened in low 
meta l licity environments ([Kudritzki et al. 19891 : IVink et al.l 
l200ll : iMevnet et al1l2006l : iMokiem et al.ll2007f ). resulting in 
lower mass and angular momentum losses from the binary 
systems in which they are present. Such angular momentum 
losses may be reflected in both the final mass of the compact 
object and the evolutionary path of the binary, i.e. whether 
or not the binary is disrupted by the supernova kick. How- 
ever, on its own, the reduced metallicity of the SMC cannot 
explain the large number of HMXB systems in this galaxy 
jDravl [20061. A recent increase in the star formation rate, 
such as required bv lHarris fc Zaritskvl l|2004 ). possibly due 
to tidal interaction between the SMC and its nearest neigh- 
bour the Large Magellanic Cloud (LMC), is necessary to ex- 
plain this exc ess of HMXBs. This i s further substantiated by 
the results of iGrimm et al.l (|2003l ), which have shown that 
the number of HMXB systems in a galaxy can be related 
to the galactic star formation rate. Using indicators such as 
the far infrared, Ha an d ultraviolet to predict the nu mber 
of HMXBs in the SMC (|Shtvkovskiv fc Gilfanovll2005l L one 
finds that star formation rate alone cannot account for the 
number of HMXBs discovered so far in the SMC. Hence, it 
seems likely that the increased star formation rate in the 
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Table 1. The spectral types of the Milky Way Be/X-ray 
binaries expanded from iNegueruelal (Il998h ■ 



07 08 09 BO B1 B2 B3 B4 B5 B6 B7 B8 B9 
Spectral Type 

Figure 1. The spectral distribution of Milky Way Be/X-ray bi- 
naries (solid) and the Milky Way isolated Be stars (dashed). 



SMC combined with the reduced metallicity environment 
have given rise to the large observed population of HMXBs. 

With the advent of arcsecond resolution X-ray tele- 
scopes the number of optically identified Be/X-ray bina- 
ries (all but one of the HMXBs in the SMC are Be/X-ray 
binaries) in the SMC has risen dramatically over the last 
few years. As there are clear differences in the numbers of 
HMXBs between the Milky Way and SMC, which can be as- 
cribed to metallicity and star formation, there may be other 
notable differences in the populations. Most fundamentally, 
how do the metallicity and star formation rate reflect on 
the spectral distribution of the optical counterparts to the 
Be/X-ray binary population of the SMC? Although their X- 
ray properties are well-studied, the spectral classifications 
of only five of the B e /X-ray binaries i n the SMC are know n 
|Covino et al l l200ll ; Icoe et al.l |2002| ; ISchurch et all 120071 ). 
An accurate spectral distribution of SMC sources can pro- 
vide constraints on the evolutionary models put forward to 
explain the spectral distribution of Be/X-ray binaries. In ad- 
dition, it may reflect the lower metallicity of the SMC and 
thus provide valuable insights into both binary evolution 
and star format i on in a metal-poor environment. 

INeeueruelal (|l998h showed that the spectral distribution 
of Be stars occurring in Be/X-ray binary systems is signifi- 
cantly different from that of isolated Be stars in the Milky 
Way. Whereas isolated Galactic Be stars show a distribution 
beginning at the early B-types and continuing through until 
AO, the Be star companions of X-ray binaries show a clear 
cutoff near spectral type B2. Figure[T]shows the spectral dis- 
tribution of Milky Way isolated Be stars (dashed histogram) 
compared with the Milky Way Be/X-ray binaries (solid his- 
togram). The Be/X-ray binary data used in the histogram 
are shown in Table [T] while the d ata used fo r the d istribu- 
tion of isolated Be stars are from Slettcbadj (ll982T l with a 
magnitude limit of V ^ 6. Such a magnitude cutoff basi- 
cally preselects the early-type Be stars and thus we expect 
that the true distribution may peak towards later-types than 
shown in Fig[T] 



2 SAMPLE SELECTION 

The objects selected for spectroscopic follow up in this work 
have been selected from a sample of SMC X-ray pulsars 
which have positional uncertainties small enough (generally 



Name 


Spectral Type 


Ref 


4U 0115+634 


B0. 2 V 


(1) 


RX J0146.9+6121 


Bl V 


(2) 


IGR J01363+6610 


Bl 


(3) 


LSI+6I 303 


B0 V 


(1) 


V0332+53 


08-9 V 


(5) 


4U 0352+30 


BO V 


(6) 


RX J0440. 9+4431 


B0 V 


(7) 


A 0535+262 


B0 III 


(8) 


SAX J0653+0533 


B2 V - Bl III 


(9) 


MXB 0656-072 


09.7 V 


(10) 


4U 0726-26 


08-9 V 


(11) 


RX J0812.4-3114 


B0.5 V 


(12) 


GS 0834-43 


BO-2 III-V 


(13) 


RX J1037.5-5647 


B0 V-III 


(14) 


A 1118-616 


09.5 V 


(15) 


4U 1145-619 


BO. 7 V 


(16) 


GX 304-1 


B2 V 


(17) 


2S 1417-624 


Bl V 


(18) 


RX J1744.7-2713 


B0.5 


(19) 


Cep X-4 


B1-B2 


(20) 



Negueruela fc Okazak] [|200lll; (21 iReig et al. 
Reig et alj J2005h (4 ) iHutchings fc Cramptoii 



iLvubimkov et al.l (119971 1 

Steele et al. I \i9m -. 



(1) 
(3) 

(5) INegueruelal [1199811 : (6) 
(7) iMotch et al] lll997ll ; (8_ 
Kaaret et al] dl999ll:~ (10) IPakull et all 
Negueruela et al] 119961) ; (12 
Israel et al] d200(Jl: ~(14) IMotch et a". 
Janot-Pacheco et al] <198ll) ; (1 6) iFeast etaT 
(171 IParkes et al ' 



1997); 



1981 



191 ISteele et al] 
199S 



200: 

Motch et aU Il99'. 

dl997h 



(1961); 

1980); ( 181 iGrindlav et al] IU9 841; 




19991 1: (201 Bonnct-Bidaud fc Mouchetl 



< 5") to allow the unambiguous identification o f their op- 
tical c ounterparts. The objects are drawn from ICoe et ahl 
(2005), with the exception of those that have no known well- 
defined X-ray position. We also exclude SXP0.72(=SMC X- 
1), which we know to be a Roche- lobe overflowing supergiant 
system, and SXP 8.02, which is now kn own to be an anoma- 
lous X-ray pulsar (|McGarrv et al.l [ 20 05). In addition to these 
objects, we include in our sample recently discovered X-ray 
pulsar systems with small positional errors, or previously 
known systems which have refined positions through detec- 
tion with either Chandra or XMM-Newton or through pre- 
vious optical follow up. Table [2] lists the sample, along with 
the method used to determine the optical counterparts and 
any comments relating to uncertainties in the identification 
of the correct optical counterpart. 

The resulting sample comprises 37 stars with V- 
magnitudes in the range 14-17. 



3 OBSERVATIONS AND DATA ANALYSIS 

With a view to characterising the spectral population of Be 
star counterparts of X-ray binaries in the SMC, we acquired 
spectra of 32 of the objects listed in Table [2] which do not 
yet have spectroscopically determined spectral types. 

The data were acquired on the nights of 2006 Septem- 
ber 13-15 and 2007 September 18 and 19 using the EFOSC2 
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Table 2. Sample of SMC Be/X-ray binary systems with 
well-defined optical counterparts. The source name in col- 
umn 1 is based on the pulse period. The column labeled 
T indicates the telescope used to determine the best posi- 
tion of the source: C - Chandra, X - XMM, A - ASCA, R 
—ROSAT, and O — refers to optical identification by pho- 
tometry or spectroscopy of objects within the error circle. 
Column 4 gives a literature reference for X-ray error circle 
or optical identification. 



. " l J \ 1 1 ^ V 


Name 




T 


Kef 


SXP0.92 1 


PSR J0045-7319 




O 


(1) 


SXP2.37 


SMC X-2 




o 


(2) 


SXP2.76 


RX J0059.2-7138 




o 


(3) 


SXP3.34 


RX J0105.1-7211 




A 


(4) 


SXP6.85 


XTE J0103-728 




X 


(5) 


SXP7.78 


SMC X-3 




C 


(6) 


SXP8.80 


RX J0051. 8-7231 




K 2 


(7) 


SXP9.13 


RX J0049.5-7310 




R 


(4) 


SXP15.3 


RX J0052.1-7319 




O 


(8) 


SXP22.1 


RX J0117.6-7330 




O 


(9) 


SXP31.0 


XTE J0111-7317 




O 


(8) 


SXP34.1 


RX J0055.4-7210 




C,0 


(6,22) 


SXP46.6 


1WGA J0053. 8-7226 


O 


(10) 


SXP59.0 


RX J0054.9-7226 




X 


(11) 


SXP65.8 


CXOU J010712.6- 


723533 


c 


(12) 


SXP74.7 


RX J0049.0-7250 




o 


(13) 


SXP82.4 


XTE J0052-725 




C 


(6) 


SXP91.1 


RX J0051. 3-7216 




K 


(14) 


SXP101 


RX J0057.3-7325 




c 


(12) 


SXP138 


CXOU J005323.8- 


722715 


C 


(6) 


SXP140 


2E 0054.4-7237 




X 


(11) 


SXP152 


CXOU J005750.3- 


720756 


c 


(15, 11) 


SXP169 


2E 0051.1-7214 




K 


(16) 


SXP172 


AX J0051. 6-7311 




X 


(17) 


SXP202 


XMMU J005920.8 


-722316 


X 


(18) 


SXP264 


AX J0047.3-7312 




X 


(17) 


SXP280 


RX J0058.8-7272 




X 


(11) 


SXP304 


RX J0101. 0-7206 




C 


(15) 


SXP323 


AX J0051-733 




X 


(17) 


SXP348 


2E 0101.5-7225 




c 


(12) 


SXP455 


RX J0101. 3-7211 




X 


(19) 


SXP504 


RX J0054.9-7245 




C 


(6) 


SXP565 


CXOU J005736.2- 


721934 


C 


(15) 


SXP700 


CXOU J010206.6- 


714115 


c 


(12) 


SXP701 


XMMU J005517.9 


-723853 


X 


(20) 


SXP756 


RX J0049.7-7323 




X 


(17) 


SXP1323 


RX J0103.6-7201 




X 


(21) 



iMurdin et all dl979l 




Coeetal 



Edge et all 



(11 Ib5J «1: 
ISouthwell fc Charles! 
(51 lHaberl fe Pietschl 
(7) lHaberl fc Sasakil <200C|); (8) ICovino et all 
(9) ICoe et all [120021 ); (10 ) iBucklev et al. " 
Sasaki et al.l 1200311; (12) iMcGowan et al 
Stevens et all l ll999l) 



(12) 
(14) 



200 
200 

Yokogawa et al. (2003 
Macomb et al.l (120031); (161 ICowlev et al.l' (Il997l 
Haberl fc Pietschl j20o3> ; (1 81 iMaiid et all [120041 
Sasaki et all [|200ll1; " (20) lHaberl et all [120041) ; 
Haberl fc Pietschl [I2005T ); (22) Coe 2007 priv. comm 




(11) 
(13) 
(15) 
(17) 
(19) 
(21) 



1-Be ctpt to radio & X-ray pulsar, 2—2" 
ambiguity in opt ctpt 



error circle - no 



faint object spectrograph mounted at the Cassegrain focus 
of the 3.6 m telescope at La Silla, Chile. We used a grism 
with lines ruled at 600 mm" 1 and which covered the wave- 
length range A A 3095-5085 A. Using a slit width of 1" and a 
CCD binning factor of 2 x 2 throughout resulted in spectra 
at a resolution of 6 A. For fainter stars we used a lower dis- 
persion grism ruled at 400 mm" 1 , covering the wavelength 
range AA3050-6100 A and resulted in a spectral resolution 
of ~ 10 A. Exposure times were adjusted for the stellar mag- 
nitude and observing conditions, with the average exposure 
time being around 1500 s. This was sufficient to achieve an 
average signal-to-noise ratio of > 100 per pixel at ~ 4000 A. 
Wavelength calibration was performed by taking compari- 
son spectra of Helium and Argon lamps through the same 
instrumental setup. Images of an evenly illuminated surface 
within the telescope dome were used for flat-fielding. 

Data reduction, comprising bias subtraction, flat- 
fielding and wavelength calibration, was undertaken with 
version 2.12.1 of the Image Reduction and Analysis Facility 
(IRAF) provided by the NOAO. With many of our exposures 
being 30 to 4 minu t es in duration, we applied an algorithm 
developed by iPvchl (|2004h to remove cosmic rays from the 
spectral images. Object spectra were then inspected by eye 
to determine the location of the target on the CCD, after 
which a low order polynomial was employed to trace the 
spectrum across the CCD. All spectra were e xtracted usin g 
the optical extraction algorithm presented bv lHornel |[l986). 
Comparison spectra were extracted along the same trace as 
the individual science spectra they were used to calibrate. 

Although flux calibration was not desired for this 
project (spectral classification was to be performed by iden- 
tifying lines rather than fitting atmospheric models), a single 
spectrophotometric standard star (LTT 7987) was observed 
to get an idea of the rough response of the telescope, grism 

and CCD as a function of wavelength. 

Finally, spectra were shifted by — 150kms - 1 (|Allenl 

Il975t) to account roughly for the recession velocity of the 
SMC and hence to place spectral features at approximately 
the correct wavelengths. 



4 CLASSIFICATION 

Although many atte mpts have been mad e to apply the 
Morgan-Keenan (MK. [Morgan et al.lll943T l sys tem of stel- 
lar classification to stars observed in the SMC l|Humphrevsl 
1 19831 ; lAzzopardi fc Vigneaulll975t ). one typically comes up 
against two problems. Firstly, the low metallicity of the SMC 
means that the metal absorption lines are very weak, making 
the classification of B-type stars, which rely on metal-helium 
ratios, difficult. Secondly the MK system was devised from 
standard stars in the Galaxy. Classification of SMC stars 
straight onto the MK grid leads to contradictions between 
the hydrogen and metal lines and translating the MK system 
through a change in metallicity is a non-trivial task. 

The defining characteristics of B-stars in the MK system 
are that they show neutral He li nes in their spectra, but n o 
ionised He (as found in O-stars Ijaschek fc Jaschekl Il990). 
with the maximum strength of the He I lines being reached 
around spectral class B2. Hydrogen line strengths peak at 
A2, so the H and He lines show opposing trends through the 
B-type sequence. The metal lines of Si II, III, IV, Mg II and 
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Table 3 . Clas si fication crite r ia for B-type stars in the SMC, from 
iLennonl (|l997l l: lEvans et al.l l|2004h . 



Criterion 



He II A4200 ~He I A4143 09 

He II AA4686, 4541 present, A4200 weak BO 

He II AA4200, 4541 absent, A4686 weak BO. 5 

He II A4686 absent, Si IV AA4088, 4116 present Bl 

Si IV A4116 absent, Si IV A4088 <0 II B1.5 

Si IV, Si II absent, Si III A4553 > Mg II A4481 B2 

Si III A4553 ~ Mg II A4481 B2.5 

Si III A4553 < Mg II A4481 B3 

Si III absent, Si II A4128/4132 < He I A4121 B5 

He I A4121 < Si IK He I A4143 B8 
Mg II A4481 < He I A4471 

Mg II A4481 > He I A4471 B9 
Fe II A4233 < Si II A4128/4132 



ratios of these to He I (and to a lesser extent the C II, O II 
and N II) are used for spectral classification of B-type stars. 

In early B-types the N II A3995/He I A4009 and 
Si IV A4089/He I AA4026,4121 line ratios, which become 
larger towards more luminous stars, can be used to distin- 
guish luminosity effects. Also, the He I A4121/He I A4144 
ratio, strengthening towards more luminous stars, is used 
for luminosity class determination. Towards later B-types 
the profiles of the hydrogen lines are used, with the line pro- 
files being narrower out of the main- sequence. It is worth 
pointing out that, for Be stars, use of the Balmer lines in 
stellar classification is limited due to the frequent infilling of 
these lines by the circumstellar disc. 

To over come th e setb acks of applying the MK system 
in the SMC. ILennonl (|l997l ). using high signal-to-noise spec- 
tra of SMC supergiants, devised a system for classification 
of stars in the SMC. This system is "normalised" to the 
MK system such that, moving through the range of spec- 
tral types, stars in both systems exhibit the same trends 
in their line strengths. This cla ssific ation method was im- 
plemented by lEvans et al. | (|2004h and lEvans et al. I |2006l ) in 
recent spectroscopic surveys of massive stars in the SMC, 
the LMC and the Galaxy. For classification of the spectra in 
thi s work, we hav e used the tempe rature criteri a as se t out 
by ILennonl (|l997T ) and utilised by lEvans et all l|2004h (see 
Table EJ. 

Wherever possible, we used the trends in the Si/He and 
H/He line ratios to estimate the luminosity classes of the 
observed stars. However, in many cases the Si lines were not 
evident due to the combined effects of low metallicity, ex- 
treme rotational line broadening intrinsic to Be stars and 
the 6 A spectral resolution of the data. Furthermore, it was 
difficult to tell whether the He lines were suffering from the 
same emission veiling as manifest in many of the Balmer 
lines. So, for the purposes of cross-checking these luminosity 
class estimates, we used the apparent V-magnitudes of the 
Be st ars together with th e distance modulus of the SMC 
(18.9, lHarries et al.ll2003T ) to determine whether the abso- 
lute magnitude of the Be star was consistent with the es- 
timated luminosity class at the given sp ectral type (using 
absolute magnitudes of OeBe stars from IWegner 2006). It- 
erating through these procedures of firstly identifying the 
spectral type, secondly using spectral features to determine 




0.8 - 



4000 4200 4400 4600 4800 5000 
Wavelength (A) 

Figure 2. The spectrum of SXP701, with defining features 
marked on the plot. 



the luminosity class and finally validating this luminosity 
class against the absolute magnitude resulted in a consis- 
tent set of spectral and luminosity classifications for these 
stars. 

This method is illustrated by the spectrum of SXP701, 
which appears in Fig. [21 The presence of ionised He in the 
spectrum, at the wavelengths A4686 A and AA4200, 4541 A 
allows us to determine the spectral type as 09.5. The lu- 
minosity criteria at this spectral type are the line ratios of 
Si IV A4089/He I AA4026, 4121, 4144 and Si II A4116/He I 
A4121. As the Si lines are weakened by the lower metallicity 
environment of the SMC, we use the strengthening He II 
A4686 alongsid e the C III/O II blend to estim ate a luminos- 
ity class of V dWalborn fe Fitzpatrick| [T990'l. The apparent 
magnitude of V = 15.87 confirms this luminosity class, as 
a star of higher luminosity class would need to be of later 
spectral type, and this can certainly be discounted on the 
basis of the He II lines. 

For fainter stars which were observed with the lower 
dispersion grism, we employed two methods that would give 
us an idea of the general spectral type of the star, rather 
than attempting to use very broad and shallow spectral lines 
to pinpoint the exact spectral type. The dereddened spectra 
(using an extincti on of E(B — V ) = 0.1(J3) were compared to 
template spectra {Pickles 1998) of different spectral types, 
using the general continuum and line features to constrain 
the spectral type. In addition, the He I AA4009, 4026/He line 
ratio, which decreases rapidly through the B main sequence 
and is hardly influenced by emission from the circumstellar 
disc, was used to gauge the rough spectral class. Once again, 
the apparent magnitude of the star was used to confirm 
the spectral and luminosity class. Figure [3] demonstrates the 
value of this approach for SXP138. It is clear just from the 
slope of the spectrum that the star is earlier than B3, but 
probably later than B0 due to the lack of He II absorption 
features in the spectrum. 



1 This is the average extinction for optic al counterparts of Be/ X- 
ray systems in the SMC as measured by IZaritskv et alj l l2002fl at 
http: / /ngala. as. arizona.edu/dennis/smcext. html 



Spectral type 1.2 
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Table 4. Properties of SMC Be/X-ray binaries. So urce names are relat ed their X-ray pulse 
peri ods. In the Ref column, refers to colours fro mlUd alski et al.Nl998h, M refers to colours 
from lMassevl l|2002h . while Z refers to colours from lZaritskv et al Tl |2002h . The Grism column 
refers to the ruling, in 1/mm of the grism used in each case. The Spec column lists the clas- 
sifications determined from blue spectra in this work unless otherwise noted. The luminosity 
class, as determined from the spectrum in the case of 600 1/mm data, and as determined 
from the absolute magnitude in the case of 400 1/mm data, of each source is given in the hum 
column. 



Source 


\ ' 


AV 


B — 


V 


A(B — V) 


Ref 


Grism 


Spec 


Lum 


cvDn no 
bA.rU.yz 


16.18 


0.02 


-0. 


.21 


0.03 


U 


600 


B0. 5— B2 


IV— V 


oArz.o i 


1 ft QO 
ID. OO 


n no 
U.UZ 


n 

U. 


n o 
.uz 


n c\a 
U.U4 


rj 

li 


ftnn 
bUU 


L^y.o 


TTT \T 
111— V 


QVDO "7 ft 


14.01 


0.08 


0. 


.06 


0.09 


V 


600 


Bl— B1.5 


TT TTT 
11—111 


QYDQ QA 
bA.ro. o4 


1 ftQ 
lo.bo 


n riQ 
U.Uo 


n 
-U 


n i 
.Ul 


n nf; 
U.UO 


f\ 

U 


ftnn 
bUU 


pi "QO 

rSl— rSZ 


TTT \T 
111— V 


bA.rb.OD 


14.59 


0.02 


-0. 


.08 


0.02 


M 


600 


<Jy.D— rSU 


IV— V 


QVD7 "7Q 

bAr 7 . Jo 


14.91 


0.02 


0. 


.00 


0.03 




600 


Bl— B1.5 


IV— V 


ovdo on 
bA.T5.oU 


14.87 


0. 12 


-0. 


.27 


0.13 


U 


600 


Uy. O—rSU 


IV— V 


bA.ry.lo 


16.51 


0.02 


0. 


.01 


0.04 


u 
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5 RESULTS 

The spectral distribution of SMC Be/X-ray binaries is shown 
in Fig. [4] with individual spectral types given in Table [4] Our 
procedure for the binning of stars where we were unable to 
constrain the spectral type to one subtype involved allo- 
cating the object fractionally to the spectral classes within 
the estimated spectral range of the source. For example, 
SXP82.4 which can be constrained to the B1-B3 range 
translates into the histogram as i Bl, 1 B2 and | B3. 

The distribution shows a similarity to t he spectral 
distributions of the Galactic l|Negueruelal Il998l ) and LMC 
<|Negueruela fc Cod [20 02 ) Be/X -ray binaries. The spectral 
distribution of Be/X-ray binary counterparts in the SMC 
peaks at spectral type Bl, compared to the LMC and Galaxy 



distributions, which peak at B0. The Galactic and LMC dis- 
tributions show a sharp cutoff at B2, whereas there are 5 
SMC objects with possible spectral types beyond B2, but 
as one can see from Table 0] the exact spectral type cannot 
be determined with certainty in these cases. A Kolmogorov- 
Smirnov test of the difference between the SMC and Galac- 
tic distributions gives a K-S statistic D — 0.22, indicating 
that the null hypothesis (which is that the two distributions 
are the same) cannot be rejected even at significances as 
low as 90%. Hence, it is likely that both Galactic and SMC 
Be/X-ray binary counterparts are drawn from the same pop- 
ulation. 

The most noteworthy selection effect, which may mask 
the real distribution of spectral types in the SMC, is magni- 
tude related: fainter stars, which are most likely to be those 
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Figure 3. The spectrum SXP138 plotted between the template 
spectrum of a BO V star and the template spectrum of a B3 V 
star. 
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Figure 4. Spectral distribution, as determined from high signal- 
to-noise blue spectra, of Be/X-ray binaries in the SMC (dot- 
dashed) as compared the distribution of Be/X-ray binaries in the 
Galaxy (solid). 

of later B-types, are the most difficult to observe and clas- 
sify. Table [4] shows that the latest spectral types observed 
among the SMC systems are those of SXP74.7 (~B3) and 
SXP101 (B3-B5). Although SXP74.7 is faint, spectral types 
of stars much fainter than SXP101 have been confidently de- 
termined. Thus, there seems very little evidence that mag- 
nitude related selection effects have influenced our sample. 
A second factor reduces the number of optical counterparts 
included in our sample in the first place: this is related to 
those SMC X-ray pulsars which do not as yet have an iden- 
tified optical counterpart. This is predominantly an effect of 
large uncertainties on the positions of pulsars as measured 
by non-imaging X-ray satellites, and we do not expect the 
absence of these objects in our sample to introduce spectral 
type selection effects. 



6 DISCUSSION 

A number of authors llRappaport Sz van den Heuve]|[l98i 
Ivan Bever fc Vanbeverenlll997l ) suggest the possibility that 
Be stars in X-ray binaries may be a product of binary evo- 
lution. During mass transfer the secondary star, or at least 



its outer layers, may be spun up (|PacketJll98il ). As it is well- 
known that the Be phenomeno n is closely associated w ith 
high stellar rotational velocities l|Porter fe Riviniusll2003h it 
is not far-fetched to suppose that such binary mass trans- 
fer may rejuvenate the secondary star and, somehow, give 
rise to the Be phenome non. Studies of Be stars in clusters 
l|McSwain fc Giesl l2005) indicate that most cases of the Be 
phenomenon may be brought about by spin-up processes 
in close binaries, rather than Be stars being born as fast 
rotators or being spun-up in the final stages of their main 
sequence lives. 

Such a model of rejuvenation may also account for the 
narrow spectral distrib ution of Be/X-ray bina ries. In popu- 
lation synthesis models iPortegies Zwartf l|l995l ) includes the 
effects of both supernova kick velocities and angular momen- 
tum loss through the L2 point. He finds that this angular 
momentum leakage is just as influential in restricting the fi- 
nal number of Be binary pulsars as the effect of a supernova 
kick and, in addition, that angular momentum loss restricts 
the spectral distribution of the Be star companions. 

IPortegies Zwartj (|l995l ) shows clearly that a higher loss 
of angular momentum restricts the mass of the Be star coun- 
terparts to > 8iW©, i.e. to spectral types earlier than B2 V 
(|Allenlll973T ) . This is caused by the fact that binary systems 
holding late-type stars, which have smaller mass ratios, tend 
to undergo a spiralling-in effec t when mass is lost through 
the 1/2 point l|Pols et al.lll99ll ). In these cases the compo- 
nents may merge or form a common envelope. Whatever 
the eventuality in such a case, it is clear that a Be neutron 
star system will not be the ou tcome. 

iNegueruela fc Coel (|2002h have shown that Be/X-ray bi- 
naries in the LMC follow the same distribution as those in 
the Milky Way, and this work (see Fig. Q shows that the 
spectral distribution of Be/X-ray binaries in the SMC is also 
consistent with that of the Milky Way systems. 

Although no major spectral survey of Be stars 
in the SMC has be en performed, preliminary results 
l|Martavan et al.l 120061 ) obtained from photometric colours 
of ~ 7700 emission line stars indicate that the spectral dis- 
tribution of isolated Be stars in the SMC is similar to that 
in the Galaxy. Hence, at this time we cannot exclude the 
possibility that isolated Be stars in the SMC may have a 
different spectral distribution from those in the Milky Way. 

Figure [5] shows the arbitrarily scaled spectral distri- 
bution of SMC Be/X-ray binaries superimposed on the 
predicted spectral dis tribution of Be/X-ray binaries (from 
IPortegies Zwartj Il995l ). As with the Milky Way and LMC 
distributions, the SMC distribution cuts off around spectral 
type B2 (~ 8Mq), indicating that there may be significant 
angular momentum losses in the binary system prior to the 
Be/X-ray binary evolutionary phase. A possible interpreta- 
tion of the fact that there is no significant metallicity de- 
pendence of the spectral distributions of Be/X-ray binaries 
is that the angular momentum is lost through mechanisms 
other than the stellar winds o f early-type com p onent s of 
these systems. For example, as IPortegies Zwartl (|l995l ) in- 
fers, the mass may be lost through the L2 point during the 
initial phase of mass transfer which serves to spin up the star 
to Be star status. It is also worth noting that there are no 
SMC Be/X-ray binaries with masses greater than ~ 22 Mo. 
This may be due to the fact that heavier systems go on to 
become supergiant X-ray binaries, which have lifetimes sig- 
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Figure 5. The absolute normalised number distribution of Be 
stars with a neutron star companion for four evolutionary scenar- 
ios, v represents the supernova kick velocities (in kms" 1 ) while 
/3 represents the angular momentum lost from the binary system. 
The solid histogram represents the spectral distribution of SMC 
Be/X-ray binaries f rom t his work (spectral type/mass calibra- 
tion from[Zombcck (2007), arbitra rily scaled along vertica l axis) 
superimposed. Original figure from [pprtcgics Zwar j <|l995h . 

nificantly shorter than the Be/X-ray systems and are conse- 
quently less abundant (there is only one in the SMC - SMC 
X-l). 
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